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Abstract 

I will briefly survey the most important results obtained so far on chaos 
in partial difi'erential equations. I will also survey progresses and make 
some comments on Navier-Stokes equations and turbulence. 
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1. Introduction 

In my view, the two most important problems in the area of chaos (tur- 
bulence) are 

1. An effective description of chaotic (turbulent) solutions beyond the 
Reynolds average, so that important characteristics of the chaotic 
(turbulent) solutions can be revealed accurately. 

2. Proofs of the existence of chaos (turbulence), so that the parameter 
regime and certain characteristics of chaos (turbulence) can be 
discovered. Control of chaos (turbulence) is a sub-problem here. 

So far, there is no substantial progress on problem 1 [30]. There are 
important progresses on problem 2 [31] which will be addressed later. A 
notable mathematical problem on Navier-Stokes equations is the global 
regularity problem which was chosen as one of the seven Clay problems. 

Extensive studies have been conducted on chaos in ordinary dif- 
ferential equations (ODE) [1], and intensive studies are focused upon 
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Lorenz equations [34]. Existence of chaos in Lorenz equations has not 
been proved by hand. Recently, there have been a new trend of com- 
puter proofs [36]. In general, for near integrable ODE with transversal 
homoclinics, existence of chaos can often be proved by hand. The same 
is true for PDE. 

Chaos in partial differential equations (PDE) has long been an open 
area. During the last decade, a standard program was established for 
proving the existence of chaos in near integrable PDE [25]. Around 
transversal homoclinics, existence of chaos can be proved by hand [20] 
[21] [22], while around non-transversal homoclinics, existence of chaos 
can be proved by hand up to nasty generic conditions [23] [24]. There 
are also attempts of computer proofs on the existence of chaos in PDE. 

In contrast to ODE, PDE has a lot of novelties. For instance, 
boundary condition plays a major role on dynamics. So far, all the 
chaos in PDE is proved under periodic boundary condition. In fact, odd 
or even constraint besides the periodic boundary condition is crucial. 
Often one can prove the existence of chaos under odd or even constraint, 
but not in general. Periodic boimdary condition permits a representa- 
tion of solutions by Fourier series. The sequence of Fourier modes in the 
Fourier series leads to a natural generalization of ODE to infinite dimen- 
sional systems. The Fomier series also makes hand calculation possible. 
Under other boundary conditions, hand calculation is often impossible. 
More importantly, the dynamics is dramatically different. For example, 
for the whole space problem (i.e. under decaying boundary conditions), 
existence of chaos has not been proved. Of course, the PDE type is fun- 
damental to dynamics. So far, the chaos proved is for mixed hyperbolic 
and parabolic semilinear systems which are near hyperbolic semilinear 
integrable systems. For semilinear systems, existence of invariant man- 
ifolds is easy to obtain. For hyperbolic quasilinear systems, existence 
of invariant manifolds is a big open problem. Here one of the interest- 
ing examples will be Euler equations. Even for the simple integrable 
derivative nonlinear Schrodinger equation (DNLS) 

iqt = qxx-i (IqI'^q)^ > (i-i) 

existence of invariant manifolds is open while its invariant subspaces 
are perfectly normal [32]. This difSculty is the major obstacle toward 
proving the existence of chaos when the DNLS is under perturbations. 

Chaos in PDE is a rich area, full of potential, with a crown goal 
of solving the problem of turbulence. Turbulence is governed by Navier- 
Stokes equations (NS). Turbulence often happens when the Reynolds 
number is large. Formally sotting the Reynolds number to infinity in the 
Navier-Stokes equations, one gets the Euler equations. Euler equations 
are a lot like hyperbolic quasilinear integrable systems [15] [17] [18]. 
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By posing periodic boundary condition to Navier-Stokes equations, one 
can study turbulence in the near Euler equations regime as a canonical 
problem. So far, there is no success in proving the existence of chaos in 
Navier-Stokes equations. 

2. Chaos in Partial Differential Equations 

The standard program that we developed for proving the existence of 
chaos in partial differential equations [25] involves multi-disciplinary sub- 
jects inchiding intcgrablc theory, dynamical system, partial differential 
equation, and functional analysis. The specific machineries are: (1). 
Darboux transformations, (2). Isospectral theory, (3). Persistence of in- 
variant manifolds and Fcnichcl fibers, (4). Mclnikov analysis and shoot- 
ing technique, (5). Shadowing technique and symbolic dynamics, (6). 
Specific Smale horseshoe construction. 

21. Transversal Homoclinics and Heteroclinics 

In a non- autonomous system (e.g. a periodic system), transvcrsality can 
often be signified by its Poincare period map. A typical example in this 
case for a success of proving the existence of chaos is the sine-Gordon 
system [22] [28]: 

Utt = (?Uxx -|- sinu -h e[—au + f sir? u] (2.2) 

which is subject to periodic boundary condition and even (or odd) con- 
straint 

u{t,x + 2-k) = u{t,x), u{t,—x) = u{t,x), or u{t,—x) = —u{t,x), (2.3) 

where u is a real- valued function of two real variables t and x, c and 
a > are parameters, ^<c<l,eisa small perturbation parameter, 
e > 0, and / = cost for example. 

Theorem 1. [22] [28] For (e,a) € (0,eo] x (ai,a2), where eo > 
and < oi < a2, under the odd constraint (2.3), there is a homoclinic 
orbit; under the even constraint (2.3), there is a heteroclinic cycle. In a 
neighborhood of the homoclinic orbit or heteroclinic cycle, there is chaos. 
That is, there is a Cantor set S of points, which is invariant under an 
iteration of the Poincare period map for some K. The action of F^ 
on S is topologically conjugate to the action of the Bernoulli shift on two 
symbols and 1. 

The proof on the existence of a homoclinic orbit or a heteroclinic 
cycle involves Darboux transformations, isospectral theory, persistence 
of invariant manifolds and Fenichel fibers, and Melnikov analysis and 
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shooting technique. The proof on the existence of chaos involves shad- 
owing technique and symbolic dynamics. 

The concept of a homoclinic orbit or a heteroclinic cycle can be 
generalized to a homoclinic tube or a heteroclinically tubular cycle [22] 
[28] by replacing orbits with invariant tubes. It turns out that the above 
theorem is still true for tubes. On the other hand, dynamics inside the 
tubes can be chaotic too. So we have a smaller scale chaos embedded 
inside a larger scale chaos. In principle, this process can continue to 
smaller and smaller scale chaos. Thus we can have a "chaos cascade" 



[22] [28]. For example, let / = f{t,e) in (2.2) where f{t,0) is periodic 
in i and ^ e T^. For example, 9 = {61,62,93), 6^ = w„i + 6^ + ei^^n, 
> 1, and i^n's are given by the ABC flow [5], 



which is chaotic for certain values of the real parameters A, B, and C. 
Let 6i = LOit + 6%, f can be, for example, of the form 



where w„'s {1 < n < A) form a quasipcriodic basis, 0^'s (1 < n < 4) are 
constant phases, and a and /3„ are real constants. 

Theorem 2. [22] [28] For (e,a) € (0,eo] x (ai,a2), where eo > 
and < ai < 02, under the odd constraint (2.3), there is a homoclinic 
tube; under the even constraint (2.3), there is a heteroclinically tubular 
cycle. In a neighborhood of the homoclinic tube or heteroclinically tubular 
cycle, there is tubular chaos. That is, there is a Cantor set S of tori, 
which is invariant under an iteration of the Poincare period map for 
some K . The action of F^ on S is topologically conjugate to the action 
of the Bernoulli shift on two symbols and 1 . 

In an autonomous system, a transversal homoclinic orbit or a hete- 
roclinic cycle is asymptotic to a limit cycle. A typical example in this case 
for a success of proving the existence of chaos is the Ginzburg-Landau 
equation in the near nonlinear Schrodinger regime: 



1?! = AsiniJs + Ccos'!?2 , 

'1)2 = B sin §1+ A cos dy, , 
■!?3 = C sin i?2 + -B cos i?i , 



(2.4) 
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(2.5) 




(2.6) 



where g is a complex- valued function of two real variables t and a;, e > 
is the perturbation parameter, is a real constant, and dx is a bounded 
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Fourier multiplier, 

K oo 

dxQ = — 5^ kqk sin kx , ioi q = qk cos kx , 

k=l k=0 

for some fixed large K. Periodic boundary condition and even constraint 
are imposed, 

q{t, X + 27r) = q{t, x) , q{t, -x) = q{t, x) . 

Theorem 3. [20] There exists a cq > 0, such that for any e G 
(0, eo); and \^\ > 5.8, there exist two transversal homoclinic orbits asymp- 
totic to the limit cycle qc = \ exp{— + 7]}. In a neighborhood of the 
homoclinic orbits there is chaos, in the sense of Theorem 1 with the 
Poincare period map replaced by a Poincare return map. 

The proof on the existence of the two transversal homoclinic orbits 
uses the same tools as above. The proof on the existence of chaos uses a 
shadowing lemma developed in [20] for infinite dimensional autonomous 
systems — a long time open problem. 

22. Non- Transversal Homoclinics and Heteroclinics 

In an autonomous system, homoclinics or heteroclinics asymptotic to a 
saddle is usually non-transversal. Some of such homoclinics or hetero- 
clinics can still induce chaos, e.g. Silnikov homoclinics. Nevertheless, a 
proof without any generic assumption on such chaos is still elusive. A 
typical example is the Ginzburg-Landau equation in the near nonlinear 
Schrodinger regime [23] : 

iqt = qxx + 2[|g|^ - uj'^]q + ie[qxx - aq + P] , (2.7) 

where q = q(t, x) is a complex-valued function of the two real variables 
t and x, u> G (1/2,1), a > and /? > are constants, and e > 
is the perturbation parameter. Periodic boundary condition and even 
constraint are imposed, 

q{t, X + 2n) = q{t, x) , q{t, —x) = q{t, x) . 

Theorem 4. [23] [24] There exists a cq > 0, such that for any 
€ e (0, eo), there exists a codimension 1 surface in the space of {a, (3,(jj) £ 
M+ X IR+ X M+ where to S (5, ^)/S, S is a finite subset, and oui < (3. For 
any (a, /?, uj) on the codimension- one surface, the Ginzburg-Landau equa- 
tion (2.7) has a pair of Silnikov homoclinic orbits asymptotic to a saddle. 
The codimension 1 surface has the approximate expression a = 1/k{u)) 
where k{u)) is given in Figure 1. Under certain generic assumptions, 
Smale horseshoes can be constructed in the neighborhood of the homo- 
clinic orbits. That is, there exists chaos under generic assumptions. 
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Figure 1: The graph of «; = k{w). 
23. Extra Comments 

Other success or partial success in proving the existence of chaos in- 
cludes the discrete Ginzburg-Landau equation in the near discrete non- 
linear Schrodingcr regime [10] [11] [12] [19], perturbed Davey-Stcwartson 
equation [14] [27], perturbed vector nonlinear Schrodinger equation [16]. 

Arnold diffusion can be proved for discrete nonlinear Schrodinger 
equation under certain Hamiltonian perturbation [29]. 



3. Navier- Stokes Equations and Turbulence 

The success in proving the existence of chaos in near integrablc partial 
differential equations naturally leads to the question of proving the ex- 
istence of chaos in Navier-Stokcs equations (NS) in the high Reynolds 
number (near Euler equations) regime. Like Lorenz system among sys- 
tems of ODE, chaos in Navier-Stokcs equations gathers the most interest. 
Unlike Lorenz system, high Reynolds number Navier-Stokes equations 
are a lot like near integrable systems, and Euler equations are a lot like 
integrable systems. 
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31. Lax Pairs of Euler Equations 

The 2D Euler equation can bo written in the vorticity form, 

dtn + {^,n} = , (3.8) 

where the bracket { , } is defined as 

{f,9} = {dj){dyg)-{dyf){d,g) , 

O is the vorticity, and ^ is the stream function given by, 

u = -dy^ , V = da;^ , 

and the relation between vorticity fl and stream function ^' is, 

n = da:V- dyu = A* . 

Theorem 5. [15] The Lax pair of the 2D Euler equation (3.8) is 
given as 

where 

and X is a complex constant, and (p is a complex-valued function. 

Starting from a Lax pair, one can try to develop a Darboxix trans- 
formation which is the key in the success in proving the existence of 
chaos in near integrable PDE. Some progress is made in establishing the 
Darboux transformation for 2D Euler equation. Consider the Lax pair 
(3.9) at A = 0, i.e. 

{n,p} = 0, (3.10) 
dtP + {<iJ,p} = , (3.11) 

where we replaced the notation ip by p. 

Theorem 6. [17] Let f = f{t, x, y) he any fixed solution to the 
system (3.10, 3.11), we define the Gauge transform Gf: 

P = Gfp=^^ [p^ - {d. In f)p] , (3.12) 

and the transforms of the potentials ^l and ^ : 

^ = + F , n = n + AF, (3.13) 
where F is subject to the constraints 

{n,AF} = 0, {AF,F} = 0. (3.14) 
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Thenp solves the system (3.10, 3.11) at {^,^). Thus (3.12) and (3.13) 
form the Darboux transformation for the 2D Euler equation (3.8) and 
its Lax pair (3.10, 3.11). 

For KdV equation and many other sohton equations, the Gauge 
transform is of the form [33] , 

P = Px- (5xln/)p . 

In general, Gauge transform does not involve potentials. For 2D Euler 
equation, a potential factor ^ is needed. Prom (3.10), one has 

The Gauge transform (3.12) can be rewritten as 

f f 

The Lax pair (3.10, 3.11) has a symmetry, i.e. it is invariant under the 

transform {t,x,y) {—t,y,x). The form of the Gauge transform (3.12) 
resulted from the inclusion of the potential factor ^ , is consistent with 
this symmetry. 

The 3D Euler equation can be written in vorticity form, 

dtO. + {u ■ V)0 - (O • V)u = , (3.15) 

where u = {tii,U2,u^) is the velocity, = f22) f^s) is the vorticity, 
V = {dx, dy, dz), n = V X u, and V • u = 0. u can be represented by O 
for example through Biot-Savart law. 

Theorem 7. [17] The Lax pair of the 3D Euler equation (3.15) is 
given as 

where 

L(j)= {n- V)^ , Aip= {u- V)(t) , 

A is a complex constant, and (f) is a complex scalar-valued function. 

Our hope is that the infinitely many conservation laws generated 
by A e C can provide a priori estimates for the global wcU-posedness 
of 3D Navier-Stokes equations, or better understanding on the global 
well-posedness. 

Notice that the equation (3.15) without the constraint f2 = V x 
u and V • w = is a compatibility condition of the Lax pair (3.16). 
In general. Lax pairs can support finite time blow up solutions. Our 
conjecture here is that adding only the first constraint = V x u, the 
equation (3.15) may have finite time blow up solutions which will be 
removed by the second constraint V • u = 0. 



Chaos in Partial Differential Equations 



9 



32. Invariant Manifolds and Their Zero- Viscosity Lim- 
its 

To begin a dynamical system study on Navier-Stokes equations, one 
needs to study the spectra of their linearizations at fixed points and 
their invariant manifolds. 

Following the notations of (3.8) and (3.15), we shall study the 
Navier-Stokes equations in the forms: 

dtfl + 0} = iy[Afl + f{x)] , (3.17) 

for 2D, and 

dtCt + {u ■ V)n - (fi • V)u = u[An + f{x)] , (3.18) 

for 3D where v is the viscosity, A is the Laplacian, and f{x) is the body 
force. Periodic boundary condition is imposed. 

Theorem 8. [26] When ly > 0, any fixed point of NS has invariant 
manifolds. The provable size of the invariant manifolds shrinks to zero 
as v ^ . 

The interesting open question is: Can the size of some invariant 
manifold be 0(1) as i/ ^ 0+ ? 

For simple examples of fixed points, detailed information on the 
spectra of the corresponding linear NS operator can be obtained [26] 
[13] [8]. In the 2D case, consider the simple shear fi, = cosy defined on 
the rectangular periodic domain [0, 27r/a] x [0,27r] where 1/2 < a < 1. 
The corresponding linear NS operator decouples into infinitely many 
sub-operators labeled by (^1,^2) <= Z^/{0}. 

Theorem 9. [26] The spectra of the 2D linear NS operator have 
the following properties labeled by (ki,k2): 

1. {aki)"^ + {k2-\-n)'^ > l,Vn G Z. When u ^ 0, there is no eigenvalue 
of non-negative real part. When v = 0, the entire spectrum is the 
continuous spectrum 

. a|fci| 
* 2 ' * 2 

2. fci = 0, ^2 = 1. The spectrum consists of the eigenvalues 

A = -lyn'^ , n e Z/{0} . 

The eigenfunctions are the Fourier modes 

Wnpe"^^ + c.c. , Vw„p e C , n € Z/{0} . 

As v , the eigenvalues are dense on the negative half of the 
real axis (—00,0]. Setting v = 0, the only eigenvalue is X = of 
infinite multiplicity with the same eigenfunctions as above. 
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3. ki = —1, ^2 = 0. (a). V > 0. For any a S (0.5,0.95), there is a 
unique z^*(a), 



v'32 - 3a6 _ I7a4 - 16a2 i /i 

< z>'*(q;) < 



4(a2 + l)(a2+4) ^ ' 2(a2 + 1) 

wftere the term under the square root on the left is positive for 
a G (0.5, 0.95), and the left term is always less than the right term. 
When V > i^*(a), there is no eigenvalue of non-negative real part. 
When V = v^{a), X = is an eigenvalue, and all the rest eigen- 
values have negative real parts. When v < V:t{a), there is a unique 
positive eigenvalue X{i') > 0, and all the rest eigenvalues have 
negative real parts. i'~^X{i') is a strictly monotonically decreasing 
function of v. When a G (0.5,0.8469), we have the estimate 



^RTTT - 16(a4l)(a2 + 4) ' + < ^^^^ 



where the term under the square root on the left is positive for 
a e (0.5,0.8469). 



V 8(a2 + 1) 16(a2 + i)(a2+4) " " ]] 8{a^ + 1) " 

In particular, as v 0^, Afi/) ~ 0(1). 

(b). V = 0. When a G (0.5,0.8469), we have only two eigenvalues 
Ao and — Aq, where Aq is positive, 



/q 2(1 a4(a2 +3) a^{l-a^) 

(a2 + 1) " 16(Q2 + i)(a2+4) < ^° ^ \j 8{a^ + i) ' 

T7ie rest of the spectrum is a continuous spectrum \—ia/2, ia/2]. 
(c). For afny fixed a e {0.5,0. 8A69), 

lim A(i') = Ao . 

1/^0+ 

Finally, when v — 0, the union of all the above pieces of continuous 
spectra is the imaginary axis iM. 

In general, for any fixed point of NS, the spectrum of the linear 
NS operator consists of only eigenvalues which lie in a parabolic region 
[26]. The above theorem establishes some properties of the zero- viscosity 
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limits of these eigenvalues for a simple shear. From numerical simula- 
tions on many examples [7], these eigenvalues often undergo fascinating 
deformations in the zero- viscosity limits, which can be classified into the 
following four categories: 

1. Persistence: These are the eigenvalues that persist and approach 
to the eigenvalues of the corresponding linear Euler operator when 
the viscosity approaches zero. (e.g. at 2D and 3D shears, and cat's 
eye.) 

2. Condensation: These arc the eigenvalues that approach and form 
a continuous spectrum for the corresponding linear Euler operator 
when the viscosity approaches zero. (e.g. at 2D and 3D shears, 
cat's eye, and ABC fiow.) 

3. Singularity: These are the eigenvalues that approach to a set that 
is not in the spectrum of the corresponding linear Euler operator 
when the viscosity approaches zero. (e.g. at 2D and 3D shears.) 

4. Addition: This is a subset of the spectrum of the linear Euler 
operator, which has no overlap with the zero viscosity limit set of 
the spectrum of the linear NS operator, (e.g. cat's eye.) 

Investigating the zero-viscosity limit is very important for studying 
the; NS dynamics near the Euler regime. By a combination of numerics 
and analysis in [7] , wo have tried to push the Alclnikov integral tool into 
the study on chaotic dynamics of NS. One interesting question is the 
conjecture that Euler equations have hetcroclinics [7]. 

Proving the existence of chaos in NS will provide a foundation for 
control of turbulence. Control of turbulence depends on the information 
that supports the existence of chaos. In fact, control can be regarded 
as an extra boby force in the NS equations. Thus, control of chaos 
(turbulence) is a sub-problem of existence of chaos (turbulence). From 
this perspective, investigating the existence of chaos in NS has great 
industrial value. 

33. Global Reguleirity of Navier-Stokes Equations 

A notable mathematical problem on 3D Navier-Stokes equations is the 
global regularity problem which was chosen as one of the seven Clay 
problems. Specifically, the fact [9] that 




being bounded only implies 
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being bounded for almost all t, is the key of the difficulty. In fact, 
Leray was able to show that the possible exceptional set of t is actually 
a compact set of measure zero. There have been a lot of more recent 
works on describing this exceptional compact set [2]. The claim that 
this possible exceptional compact set is actually empty, will imply the 
global regularity and the solution of the problem. The hope for such a 
claim seems slim. There are weaker conditions implying global regularity, 
e.g. L'^ norm of velocity being bounded for t in the exceptional set [4] 
[37]. But I think the weak and the strong conditions will be true or 
false simultaneously. There are several alternative ways of approaching 
the global regularity problem [35] including the Lax pair for 3D Euler 
equations in a previous subsection. 

Even for ordinary differential equations, often one can not prove 
their global well-posedness, but their solutions on computers look per- 
fectly globally regular and sometimes chaotic. Chaos and global regular- 
ity are compatible. The hallmark of chaotic solutions is their sensitive 
dependence on initial conditions. That is, a small change in the initial 
conditions will lead to a huge change after sufficiently long time, while 
all the solutions are kept in a bounded region in the phase space. This 
leads to the unpredictability, not irregularity, of chaotic solutions. The 
fact that fluid experimentalists quickly discovered shocks in compressible 
fluids and never found any finite time blow up in incompressible fluids, 
indicates that there might be no finite time blow up in 3D Navier-Stokes 
equations (even Euler equations). On the other hand, the solutions of 
3D Navier-Stokes equations can definitely be turbulent. 

Replacing the viscous term v Au by higher order derivatives, one 
can prove the global regularity [6]. This leaves the global regularity of 
a more challenging and interesting mathematical problem. Assume that 
the unthinkable event happens, that is, someone proves the existence of 
a meaningful finite time blow up in 3D NS, then fluid experimentalists 
need to identify such a finite time blow up in the experiments. If they fail, 
then the choice will be whether or not to replace the viscous term u Au in 
the Navier-Stokes equations by higher order derivatives to better model 
the fluid motion. Such a freedom was allowed in the original derivation 
of NS. Specifically, the fiuid shear stress can depend on second or higher 
derivatives of velocity besides the fluid strain given by flrst derivatives 
of velocity. 

The global regularity and turbulence are separate problems. Tur- 
bulence is more of a dynamical system problem. Of course, global at- 
tractors depend on global well-posedness [3]. Often a detailed local dy- 
namical system study does not depend on global well-posedness. Local 
well-posedness is often enough. In fact, this is the case in my proof on 
the existence of chaos in partial differential equations [25] . 
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34. Searching for an Effective Description of Chaotic 
(Turbulent) Solutions Beyond the Reynolds Av- 
erage 

Of the greatest industrial value is an effective description of chaotic (tur- 
bulent) solutions beyond the Reynolds average. From what wc learn 
about chaos in partial differential equations [25] , turbulent solutions not 
only have sensitive dependences on initial conditions, but also are densely 
packed inside a domain in the phase space. They arc far away from the 
feature of fluctuations around a mean. In fact, they wander around in 
a fat domain rather than a thin domain in the phase space. Therefore, 
averaging makes no sense at all. One has to seek other descriptions. In 
[30], an interesting "segment description" is proposed. 
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